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Abstract. We calculate the cosmic microwave background temperature bispectrum from 
cosmic strings, for the first time including the contributions from the last scattering sur¬ 
face, using a well-established Gaussian model for the string energy-momentum correlation 
functions, and a simplified model for the cosmic fluid. We check our approximation for the 
integrated Sachs-Wolfe (ISW) contribution against the bispectrum obtained from the full sky 
map of the cosmic string ISW signal used by the Planck team, obtaining good agreement. 
We validate our model for the last scattering surface contribution by comparing the predicted 
temperature power spectrum with that obtained from a full Boltzmann code treatment ap¬ 
plied to the Unconnected Segment Model of a string network. We find that including the 
last scattering contribution has only a small impact on the upper limit on the string tension 
resulting from the bispectrum at Planck resolutions, and argue that the bispectrum is unlikely 
to be competitive with the power spectrum at any resolution. 


Contents 


1 Introduction 1 

2 Temperature Fluctuations 3 

2.1 Gaussian Model 4 

2.2 Gott-Kaiser-Stebbins (ISW) effect 5 

2.3 Gontributions from Last Scattering 5 

3 Power Spectrum 7 

3.1 ISW Power Spectrum 7 

3.2 Intrinsic and Doppler Power Spectra 8 

3.3 Gomparison to GMBAGT 9 

4 Bispectrum 10 

4.1 ISW Bispectrum 11 

4.2 Intrinsic and Doppler Bispectra 12 

5 Results and Discussion 13 

6 Conclusions 17 

A Separable Decomposition 18 


1 Introduction 

Gosmic strings are linear topological defects which may be formed at symmetry-breaking 
phase transitions in the early Universe [1]. Observation of the power spectrum of the cosmic 
microwave background (GMB) radiation has gradually revealed that cosmic strings cannot 
be the dominant source of the seeding of large-scale structure in the Universe. All-sky GMB 
satellite experiments such as WMAP and Planck have shown that the cosmic strings con¬ 
tribute less than 3% of the power observed [2, 3]. Nevertheless, interest in strings remains 
very well motivated, as cosmic strings are a by-product of many particle physics models with 
extra gauge symmetries [4], and cosmic strings are generically formed at the end of brane 
inflation [5]. A detection of cosmic strings would provide important clues about the theory 
of the very early universe and, perhaps, as to the validity of the string theory paradigm. 

There are numerous methods by which the existence of cosmic strings may be tested, 
(for a review see, for example, [6-8]). However, constraints from many of these methods are 
either weaker than those from the GMB power spectrum, or subject to uncertainties in our 
understanding of cosmic string evolution. In this paper we assess the widely-held belief that 
the non-gaussianity of the cosmic string GMB signal can provide stronger constraints than the 
power spectrum. For example, in [9, 10], an edge detection (Ganny) algorithm was suggested 
as a method to search for cosmic strings in GMB anisotropy maps, suggesting that such a 
probe may prove superior to existing measures by up to an order of magnitude. 

This program has been undertaken by the authors in previous papers, specifically [11- 
14] in the context of the integrated Sachs-Wolf (ISW) signature of the GMB, and in [15] 
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for the case of matter perturbations. The latter paper established that matter polyspectra 
are unlikely to provide competitive constraints to those obtained using the CMB, while the 
former papers established that the CMB ISW bispectrum is suppressed, ultimately by the 
approximate time reversal invariance of the string network.^ 

This symmetry suppression ensures that the bispectrum of the ISW signature will not 
be competitive with power spectrum constraints. These analytic considerations were further 
established by simulations of the post-recombination signal using stacked maps [3]. The 
four point correlation function (trispectrum) of the ISW effect was also considered in [16], 
and shown not to be symmetry-suppressed. However, the relatively low signal to noise and 
difficulty of measurement means that the trispectrum is also a sub-optimal measure. 

The question remains whether the inclusion of effects from the recombination surface 
may alter these conclusions - particularly in the case of the bispectrum. It is known for the 
power spectrum that recombination effects dominates over the ISW signal for angular scales in 
the range 200 < (. < 500. The recombination surface CMB fluctuations are generated by fluid 
perturbations, and the simple arguments leading to the suppression of the ISW bispectrum 
do not apply. Hence it is important to quantify the contribution. In addition the most 
recent constraints obtained using Planck CMB data ignored this effect. Therefore, a sizeable 
recombination surface contribution may change constraints dramatically. 

Cosmic string constraints are generally characterised in terms of the dimensionless mass 
per unit length, G^, where G is Newton’s constant, and ^ represents the mass per unit length, 
which is proportional to the square of the symmetry breaking scale. Current power spectrum 
constraints on cosmic strings result in the 95% confidence bound Gjj. < (1.5 - 3.2) x 10“^ 
(with the range due to uncertainties in modelling between Abelian-Higgs and Nambu-Goto 
simulations), with the corresponding ISW bispectrum constraint giving G/r < 8.8 x 10“'^. 

The contribution to the bispectrum from the last scattering surface has thus far been 
largely ignored due to its computational complexity. While numerical approaches to com¬ 
pute the power spectrum may avail of orthogonalisation techniques (such as the Cholesky 
decomposition) to decompose the unequal time correlators (UETCs) source term into a sum 
of coherent sources, this is not feasible for the bispectrum, for which 3-point correlation func¬ 
tions of the string energy-momentum tensor would need to be computed. In order to make 
progress, we make a number of simplifying assumptions. We model the 3-point UETCs using 
a Gaussian approximation for the string network. We ignore the baryon content of the cosmic 
fluid and treat the Silk damping using a damping envelope rather than computing with a full 
Boltzmann code. Compensation effects induced by the cosmic strings on super horizon scales 
are fitted using a fitting function determined in [17] and discussed in some detail in [15, 18]. 

With our approximations, we may find the ffuid perturbations with a Green’s function 
method applied to the source functions (the model UETCs). The Green’s functions are calcu¬ 
lated numerically and the resulting (last-scattering surface) bispectrum may then be readily 
computed in an efficient manner. In order to obtain the CMB bispectrum from this bispec¬ 
trum a weighted integral over spherical bessel functions must be performed. This integration 
is practically intractable unless the bispectrum is separable, i.e. unless the bispectrum can 
be written as a sum of terms of the form f{ki)g{k 2 )h{k^). Therefore, we employ a modal 
separation technique [19, 20] which allows one to take a grid of bispectrum values and write 
it in terms of a sum of one dimensional partial waves. In a further approximation, we neglect 
correlations between the ISW signal and that from the last scattering surface, and add them 

^B. Wandelt, private communication. 
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to obtain an estimate of the fnll CMB power spectrnm and bispectrnm. In order to test the 
approximations, a comparison is made between the power spectrnm and that compnted nsing 
the fnll Boltzmann code nsing CMBACT[21]. 

Onr compntation of the ISW signatnre follows the approach ontlined in [12, 14], and is 
carried ont nsing a flat-sky approximation, with the fnll sky CMB power spectrnm and bispec¬ 
trnm recovered nsing the Limber approximation. We expect that the flat-sky approximation 
gives reasonable resnlts at angnlar scales of order 1 degree or less. As with the bispectrnm 
dne to last scattering snrface effects, the ISW bispectrnm is again decomposed into a snm 
over partial waves, allowing the calcnlation of the signal to noise to be readily compnted. 

Combination of the bispectra from recombination to that dne solely to the ISW effect 
are shown to resnlt in only a moderate decrease of the Icj error bar from 0.15 

to A(G/r/10“®)^ ~ 0.11, conhrming that indications from the ISW only calcnlation that the 
bispectrnm is a snb-optimal probe (compared to the power spectrnm) remain trne when all 
contribntions are taken into acconnt. 

The paper is ontlined as follows: In Section 2 we describe the methods nsed to compnte 
the temperatnre flnctnations dne to cosmic strings. We hrstly describe the Ganssian model 
for the compntation of the nneqnal time correlators. We then recapitnlate the calcnlation of 
the ISW signal, and ontline the method nsed in this paper to compnte the contribntions from 
the last scattering snrface. In Section 3 we describe onr compntation of the power spectrnm, 
describing the ISW spectrnm hrst, followed by the contribntions from recombination. We 
compare the total power spectrnm to that obtained nsing a fnll Boltzmann code in order to 
test the validity of the approximations nsed. Section 4 is dedicated to the calcnlation of the 
bispectrnm from both the ISW and recombination contribntions, while in Section 5 we present 
and discnss the resnlts of these compntations. Finally in Section 6 we present onr conclnsions. 
In addition we inclnde in Appendix A a discnssion of the separable basis techniqne nsed to 
project the recombination bispectra to their CMB connterparts. 


2 Temperature Fluctuations 


In the synchronons gange, given by the metric ds^ = a^(r/)[-d?]^ + (Sij + hij)dx'‘dx^], nsing 
the thin scattering snrface approximation (whereby baryons and radiation act as a tightly 
conpled flnid nntil recombination which occnrs instantaneonsly), the temperatnre anisotropy 
in direction h on the sky is given by [22] 


AT 

~T~ 


(h) = 


^(hdec) + n- v^(r/dec) - J f 
4 Z Jrt , 


VO 


drihij{ri)n^h^ , 


Vde 
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where r/dec and ryo represent the conformal times at recombination and today, respectively. 
The hrst two terms represent the intrinsic and doppler terms at the last scattering snrface 
respectively, with 6r and denoting the energy density ond velocity of the photon-baryon 
hnid. The third term, dnbbed the integrated Sachs Wolfe (ISW) effect, describes the integral 
of the metric hnctnation along the line of sight (x(ry) = x(r/dec) + n(h~hdec))) with the overdot 
denoting differentiation with respect to conformal time, ry. The sonrces of the hnctnations 
are described by their stress-energy tensor, given, in the case of a cosmic string with 

spacetime trajectory = (ry,X(fT,ry)) (where {cr,r]) denote the world-sheet coordinates), by 


0^^(x,ry) = /r J da {eX^^X^ - 


(2.2) 


where prime denotes derivatives with respect to cr, e = yX'^/(1 - X^) and we impose X-X' = 

0 . 


2.1 Gaussian Model 

Computation of the power spectrum and bispectrum of the temperature fluctuations necessi¬ 
tates the calculation of the unequal time correlators of the energy-momentum tensor. In the 
case of the power spectrum, numerical codes such as CMBACT are based on the unconnected 
segment model (USM), where the string network is modelled as a set of straight segments 
which are randomly positioned and oriented, and which move with fixed velocity at each 
conformal time, v{r]), at right angles to the orientation. Despite recent analytical advances, 
[23], in the case of the power spectrum, the calculation of unequal time three point correlators 
in the case of the bispectrum is prohibitive with this approach. Instead we use a Gaussian 
model, first described in [11], where the string network is described - for modes well inside 
the horizon - as an ensemble of randomly placed strings with a Gaussian distribution for the 
random fields X, Xb 

We denote the relevant correlation functions as: 


r(o--,r/) = U{a-,r]) = {{X{a,r]) ,r])) ■ ,r])), 

= (X((T,? 7 )-X(cr',? 7 )), (2.3) 

where we write a- = a - a'. The asymptotic small scale limit for these functions (dropping 
the explicit reference to the time dependence) is given by 

r(cj) Ri n((T) Ri Co(T/(2^), V{(7) Ri (2-4) 


with ^ denoting the correlation length, which in the scaling regime (satisfied in the deep 
radiation and deep matter eras) obeys the relation ^ ^ r]. We shall, therefore, often express 
the correlation length in the form ^(?y) = a{7])rj. The constant cq is small due to approximate 
time-reversal symmetry, and may be expressed in the form [14, 24] 


co = ^(XbW") 


2^2 1 - 81;' 
vr 1 + 81;' 


v{l - v^). 


(2.5) 


In addition the constraint X^ + X'^ = 1 implies v‘^ += 1. The velocity-one scale model is 
used to evolve these quantities subject to the equations: 


where H = dlna/dry, k = (2V2/vr)(l - 81;®)/(1 + 81;®), and c denotes a constant energy loss 
rate (with value 0.23 to agree with GMBAGT). 

While the correlation functions described in equation (2.3) are equal time correlators, 
we shall also require the unequal time correlator ([X(cr, r;) - X(cr^,r/')]^). For small time 
differences and small spatial separations (such that \ri - r]'\ = \r]^\ < ^ and |it_| < ^), we may 
approximate this quantity as [14] 

([X(a,r?)-X(^7',r?')f)-r(^7_,r/+) + D(0,r/+)r/^ (2.7) 


where r]+ = (r] + r]')l2. 
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2.2 Gott-Kaiser-Stebbins (ISW) effect 

As derived in [11], the ISW contribution resulting from a photon moving with 4-momentum 
and intersecting the past-light cone of the string, is given by 

V^—= -87rG/r J - X), (2-8) 

where u = X - (X’.pjX.p)^', and where the quantities expressed are evaluated at the light 
cone crossing time tJj. = r] + z-p-X.(a,r]r), and (5^^^ represents the two dimensional Dirac delta 
function. Using the notation 5 = 5T /T we write the expression in Fourier space in the form, 

k‘^6k{pr) = -iSnGpk^ J dau^{a, , (2.9) 

where A = 1,2 runs over the transverse coordinates. To compute the total contribution to the 
ISW effect it is necessary to integrate over all contributions from the last-scattering surface 
to today. 

The power spectrum, dP, and bispectrum, dB, contributions due to strings crossing the 
light cone at rj^ are given respectively by 

(4i(?/r)(^k2(?/r)) = (27r)^(5^^^(ki-k2)(iP(/ci,r/r), (2.10) 

(4i(?/r-)42(?/r)43(?/r)) = (27r)^(5^^^(ki + k2 + k3)dP(A:i,/i;2, fcs; 7/,.), (2.11) 

/o\ 

where is the two dimensional Dirac delta function (which imposes the triangle condition 
in the flat sky limit in the case of the bispectrum). Since the formulae expressed are already 
expressed in the flat sky one may relate the power spectrum and bispectrum to their CMB 
equivalents using the Limber approximation [25], such that we expect the ISW expression 
to be accurate for angular scales i > 100. More explicitly, the angular power spectrum and 
angular bispectrum contributions due to the string with light-cone crossing time pr in the 
form 


Ceipr) = dP{k*,pr)Ixl , 
where Xr = po - Vr: and k* = iiIXr- 


= dB{kl,kl,kl)lxl, 


( 2 . 12 ) 


2.3 Contributions from Last Scattering 

Decomposing the temperature fluctuation in terms of spherical harmonics, ^(n) = Y,im, nzmk;m(n), 
we may express the coefficients aim due to the contributions from last scattering, (2.1), as 




/ 


dn4m(n) J 


d^k r(5^(k) 


(27r)^ L 4 


+ n • V. 


.(k)]f 


zk.fiXde 


(2.13) 


where Xdec = Vo ~ Vdec- The velocity term may be written in terms of its scalar and vector 
decomposition as 


Vr{k) = iknf (k) + (k), 


(2.14) 


where .k = 0. As we shall note later, the vector part may be neglected. Therefore, the 
spherical decomposition of the temperature fluctuation may be further expanded to give 




■'/ 


d^k r6r{k) 


(27r)3L 4 


ji{kXdec) + Vr {k)j'^{kXdec) T£m(k) 
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(2.15) 









where j'^{x) = dji{x)ldx, with denoting the spherical bessel function. 

The governing eqnations we shall nse assnme a tightly-conpled photon-baryon flnid to¬ 
gether with a cold dark matter component^ [26]. The scale factor, 0 ( 77 ), evolves according to 
the Friedmann eqnation, 


da _ 




(2.16) 


We normalise the scale factor to satisfy a^rjo) = 1 today. As described in [26] this re- 
snlts in the valnes at matter-radiation eqnality for the conformal time and scale factor 
r/eq = 16.31/(nm^^)Mpc, Oeq = 1/(232190^^^)) with the conformal time today given by 
r]o Ri 13824Mpc. In addition we assnme that the conformal time at last scattering (as dehned 
nsing the peak of the visibility fnnction) is ~ 277Mpc. We nse the Planck [27] parameters, 
klm = 0.32, h - 0.67, Qa = 0.68, with = ^mOeq- 

The linear pertnrbation eqnations governing the evolntion of radiation (r) and cold dark 
matter (c) np nntil recombination are given by 


a 


^ 3/d\^ / a(5c + (2 + i2)aeqhr 

Vo/ 1 a + Oeq + ^a'^ 


47rG0+ , 


■ 4 4 ■ 

+ g V • V,. ~ 3 ~ ^ ’ 

d R 1 e 

Vr H-Vr H- VOr = 0 , 

al + R 4 + 47? 


(2.17) 

(2.18) 


where we remind the reader that overdot represents differentiation with respect to conformal 
time, and where 0+ = 0oo + 0ii is the sonrce term dne to the string energy-momentnm tensor. 
The variable, 7?, is dehned as the ratio of the pertnrbation of the baryon (5) density to the 
photon density and is given by 7? e = ^pbI^P'^ ~ 0.6(Qbh^/0.02)(a/10~^). Since the 

vector component of v^, i.e. is not sonrced, we neglect this component and consider only 
the scalar part, v^. These eqnations are solved nsing a Green’s fnnction techniqne, with 


Sc + -Sc - - I - 

a 


-m 


QjQc (2 H" 


a + aeq + 


= 47rG0+ , 


k - 4 • 

Sr - -Svr - ^Sc = 0, 

^ d R ^ n n 

^ a 1 + 7?^’'’' ^ 1 + 7?^'' " ’ 


(2.19) 


snbject to the initial conditions Q - Sr = Svr^Sc = 1 = ‘iSrj^ at p = p' and Qi\f = 0 for 
p <p' (for N e {r, s}). We have normalised the Green’s functions snch that the solntions for 
(5r./4 and vf are given, respectively, by 

^{k,p,pi) = -KG f Sr{k]p,p')e+{k,p')dp', v^{k,p,pi) = ttG f Svr{k;p,p')e+{k,p')dp'. 
4 Jr)i Jrji 


Compensation factor — One particnlar issne with this treatment is that the compensating 
nnder-density created dne to the formation of the topological defects is not inclnded. While 

■^We shall neglect the effect of neutrinos in our analysis 
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a more involved treatment may account for such effects exactly, a simple prescription was 
described in [28] is the use of a compensating factor, 7 c given by the kernel 



( 2 . 20 ) 


where kc{T])'r] interpolates between \/6 in the radiation era and \/l8 in the matter era. 


Damping factors — Our treatment of recombination has thus far assumed an instantaneous 
decoupling. In reality the last scattering surface has a finite thickness resulting in a diffusion 
damping of photons. One may account for these using a damping envelope [29], D^{k) ^ 
where for the parameters used in this paper we may set k^ = 0.109 Mpc“^ and 

= 1.45. 

We may summarise these effects in the intrinsic and Doppler terms by writing. 



We simplify notation and set 


g\k]rj,r]') = gr{k;r],r]')'yc{k,r]')V^{k) 
g'^(t,r],r}') = gr{k;r],r]’)-/cik,r]’)V^{k) 


( 2 . 22 ) 


with the superscripts i and d used to denote intrinsic and Doppler, respectively. The advantage 
of the approach used, as described earlier, is that the calculation of the bispectrum becomes 
tractable. At the level of accuracy required for this paper, the approximations made are 
deemed to be adequate to assess the relative importance of the terms at last scattering to the 
total CMB bispectrum. In §3 we shall provide a check of the qualitative and quantitative fit 
to the power spectrum using these approximations by making a comparison to that obtained 
using a full Boltzmann code treatment. 

3 Power Spectrum 

In this section we compute the CMB power spectrum due to a network of cosmic strings. 
We shall first present our computation of the ISW power spectrum, and, subsequently, our 
calculation of the power spectra due to the intrinsic and Doppler contributions from the last 
scattering surface. Due to the approximations made throughout this work, we expect that the 
power spectrum deviates on superhorizon scales on the last scattering surface. Nevertheless, 
our aim is to accurately model effects on scales of the correlation length of the string network 
(at last scattering) and below, since the signal to noise is dominated on those scales. We shall 
test the accuracy of our approximations by comparison to the power spectrum computed 
using CMBACT. 

3.1 ISW Power Spectrum 

In order to compute the total ISW power spectrum we integrate up all contributions from 
recombination to today. We shall assume that the contributions from strings at unequal 
times are uncorrelated. In the flat-sky approximation used to compute the ISW temperature 



fluctuations, the power spectrum contribution, dP, due to strings crossing the light cone at 
r] is given by equation (2.10). The total angular power spectrum, Ci is then given by (c.f. 
equation (2.12)) 

fCi= r° , (3.1) 

JrjatBxt dr] 

where k* = i/{r]Q - r]). On sufficiently small angular scales ?7start is given by r/deo while more 
generally we shall dehne rjsturt to ensure validity of the Gaussian model. More precisely the 
Gaussian model is assumed to be valid for length scales, k~^, satisfying k~^ ^ ^ = oti]. Therefore 
we dehne, T^startj such that 

r/start = max(7/dec, (^a)"^) • (3.2) 


Our calculation of the power spectrum requires the calculation of the two point functions 
of the string network functions, as described in [11, 12, 14] resulting in the power spectrum 
contribution due to strings crossing at r] given by 


dP{r]) = 




(3.3) 


where {u^{a)u^{a')) = {6ABl‘^)U{a - a'), k"^ = k"^k^, A = (27r)^h^^(0) is a formal area 
factor, and P represents the total string length of the network at t]. Since the total length 
per unit volume of the string network is given by 1/^^, we infer that 


dC c 
Adrj ar]'^ ’ 


(3.4) 


where c is a constant which was deduced to be approximately unity [14]. We use the asymp¬ 
totic small scale limit for the two point function U{a) ~ u^, which is related to v (given by 
equation (2.4)) via 


U(a) -v^ + -{1 - + v^). 
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The total angular power spectrum is therefore given by 

1 \/^U^ 1 


fCi = {8TTG]if dr] 

"J ??start 


^7start 'f] k 


erf 


k*tar] 


(3.5) 


(3.6) 


where we simplify notation by suppressing the implicit time dependence of the quantities 
k* = il(r]o - r]) and of the string network parameters u,i,v,ct- 


3.2 Intrinsic and Doppler Power Spectra 

The Fourier space power spectrum of a three dimensional quantity, h(k), is given by 

{6{ki)S{k2)) = (27r)3<5g^(ki + k2 + k3)P{k ), (3.7) 

( 3 ) 

where represents the three dimensional Dirac delta function. We shall denote by P* 
the power spectrum of the intrinsic perturbation (5^/4, and by P'^ the power spectrum of 
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the Doppler perturbation . In addition we shall denote the cross spectrum by (with 
ptxi ^ pt pdxd - pdy Following the Green’s function technique described in §2.3, the 
respective power spectra are given by 


T^axbn\ r ^^2 J ^(0+(k,??l)0+(k,?72)) 

P (fc) = (ttG) // drjidr] 2 y {k-,r]dec,'ni)y {k-,'ndec,'n 2 ) --, (3.8) 

JJrji V 

where a, 6 e {i,d}, and V = (27r)^5^^(0) represents a formal volume factor, and r]i denotes 
the initial time (taken in this paper to be 0.01 Mpc). The unequal time correlator (UETC), 
(0+0i()/V, is given in the Gaussian model by [15] 


(0+(k,7?i)0:((k,7?2)) 

V 


= V 


71 
6 a^t krjiri2 



k ‘^ v ‘^ ril 2 2(7/1772)^'^^ 

6 ) ril + ril ' 


where 7712 = - ? 72 - The final factor was advocated in [15] to account for the super-horizon 

behaviour of the UETG. 

The angular power spectrum, is defined by 



(3.9) 

and may be inferred from equation (2.15) to be of the form 


ch-j 

f P\k){j,{kx^,,)fk^dk, 

(3.10) 


f P'^{k){j[{kxAec)fk‘^dk, 

(3.11) 

^ixd _ 2 

~ j 
TT J 

f P'-'''^{k)ji{kXdec)jeikXdec)k‘^dk . 

(3.12) 


3.3 Comparison to CMBACT 

The various approximations made have been done so since the scope of this paper is to provide 
an estimate for the relative magnitude of the bispectrum due to last scattering surface effects 
compared to that due to the integrated Sachs Wolfe effect. We compute the ISW, intrinsic 
and Doppler terms using the Planck parameters as given below equation (2.16) and using the 
initial conditions for the velocity-one scale model as used by GMBAGT, i.e. v{r]i) = 0.65, 
Oi{rji) = 0.15, with rji = 0.01 Mpc. We have set the string tension to G 77 = 10^® for reference. 
In order to assess the impact of the various approximations we compare our computation to 
that obtained using a full Boltzmann treatment, including the effects of baryons and photon 
diffusion exactly. GMBAGT uses such a treatment, with the UETGs of the string network 
obtained using simulations of the unconnected segment model. While the cross term between 
the intrinsic and Doppler terms is found to be negligible, it is found that on length scales 
around the correlation length at last scattering (corresponding to angular scales I ~ 400) the 
contributions due to last scattering surface effects exceed those due to the total ISW effect. 
In Figure 1 we plot the comparison between GMBAGT, the total obtained using the Gaussian 
model approach used in this paper (including ISW, intrinsic and Doppler effects), and also, 
for reference, plot the ISW only contribution. We observe that on scales for angular scales 
i > 200 the fit to the power spectrum obtained using the approach adopted in this paper 
replicates both qualitatively and quantitatively the features of the power spectrum obtained 
using the full Boltzmann code. The discrepancy for scales i < 200 is possibly due to the 


9 










Figure 1: Comparison between the angular power spectrum obtained using the full Boltz¬ 
mann code (using CMBACT) and that obtained using the Gaussian model approach used in 
this paper. The quantity is plotted in units of the string tension, G^. We plot the total 
contribution including the integrated Sachs Wolfe, intrinsic and Doppler effects, as well as the 
contribution due to the ISW term alone. The Gaussian model is seen to provide a reasonably 
accurate fit on angular scales £ > 200. 


neglect of super-horizon physics and compensation effects in the integrated Sachs Wolfe term, 
which - calculated using a flat-sky formalism - could only be expected to be accurate for 
angular scales well above t ~ 100. We conclude that our approach is precise to a sufficient 
level of accuracy to faithfully estimate the signal to noise improvement to the bispectrum 
that might be obtained via the full inclusion of terms from the last scattering surface. 

4 Bispectrum 

The GMB bispectrum is given by the three point function of the the spherical harmonic 
coefficients of the temperature ffuctuations, 

where the G"', representing the Gaunt integral, imposes the triangle condition in multipole 
space and is given by 

= J c^ny£imi(A)y£2m2(n)F4m3(fi). (4.2) 
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One is primarily interested in the signal to noise measure in order to assess the detectability 
or otherwise of this signal, given by 



/sky ^ 


6 


'^£1 


(4.3) 


where = 

Wigner 3-j symbol (which, inherited from the Gaunt integral, imposes the triangle condition), 
and /sky denotes the fractional sky coverage of the survey considered. The power spectrum 
is the CMB power spectrum, observed to show high consistency with the inflationary 
paradigm. In order to compute the signal to noise accurately we include the effect of beam 
and noise effects in our estimates, as described in §5. Computing the signal to noise for a 
reference value of we note that Icr error bars on the quantity (G/r/10'®)^ is given by 

A(Gf./10-«)“ = , / . (4.4) 

\J{S/N)g^=io-6 

In this section we shall outline how the CMB bispectrum of the ISW, intrinsic and Doppler 
terms are computed in the Gaussian model framework described in §2. 


(2G + 1)(2£2 + 1)(24 + 1) 

dvr 


/ h h h ^ 

I 0 0 0 I’ 


with the second factor denoting the 


4.1 ISW Bispectrum 

As in the case of the power spectrum, we compute the ISW bispectrum by integrating up all 
contributions between last scattering and today, assuming that contributions from unequal 
times are uncorrelated. The computation is carried out using the flat-sky approximation, with 
the bispectrum contribution, dB due to strings crossing the light cone at rj given by equa¬ 
tion (2.11). The Limber approximation, (2.12), may then used to give the CMB bispectrum 
result. Explicitly we compute. 


rm 

'^^?start 


1 

{rjo - 


dB{k \, , ^3; r/) 


dr] 


(4.5) 


where we shall define r/gtart later in the section to ensure validity of the Gaussian model for 
the length scales integrated over, k* = iil{r]Q - rj). In the Gaussian model we may readily 
compute the contribution, dB, in a similar fashion to the power spectrum as 


dB{ki,k2,k3;r]) 


. (SttG/u)^ 
^ Akfk'^k'^ 


y" daida2das{ki k2 k!^u^{cji)u^{a2)vP , 


(4.6) 


where we set fcs = |ki +k 2 |, A = (27r)^<5^^(0), A,B,C e {1,2} and D = X]f=iki.Xj. Utilising 
the Gaussian model functions described in §2.1 this equation results in the expression [12, 14] 


dB{ki,k2,ks;r]) 

dr] 


(SttG^)" 


CqU^ 


vp{kik2k^Y 


J dai2dcri3 [klK23(^h + + fc|Ki2i732] 


xexp — 


Zi<j PKijaf. 


(4.7) 


where aij = ai -ctj, and Kij = -kj.kj. In [12] it was shown that this integral can be computed 
analytically in a sufficiently small angle limit, which may be deduced from [14] to be given 
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by the condition ~ noting that the area of the triangle defined by 

wavennmbers ki,k 2 ,k 3 is given by ■^k\k\ - While one may prodnce analytic approxi¬ 

mations to the bispectrnm on smaller scales, the signal is snppressed, so we instead consider 
only time scales satisfying this condition, i.e. we impose the start time of onr integration to 
be given by 


T^start = max 



4: kj ^ 

ai^k‘ik‘l-KiJ ■ 


(4.8) 


In £ space this translates to the condition r/start ~ max(7ydec) %(30/.^i)), where ii = £ 1^21 sin 
where 612 is the angle between the triangle sides of length ki and k 2 - Thns for triangles that 
are too small or angnlar scales which are too small snch that ii there is no contribntion to 
the bispectrnm integral. This is to be expected, since onr nse of the flat-sky approximation 
is not anticipated to be valid for angnlar scales i < 100. 

The bispectrnm contribntion is then given by 

dB{ki,k2,k3;r]) _ dTrcpif + ^ 1^12 

dr] rf‘{kik 2 k 2 ,Y (ki2«^13 + ^ 12^23 + ki3K23)^/^ 

One may readily compnte the CMB ISW bispectrnm from cosmic strings for particnlar con- 
fignrations nsing the above expressions. From eqnations (4.8) and (4.9) one may dednce that 
in the eqnilateral confignration for which ?7start ~ 30(14000/£), the ISW bispectrnm will scale 
as for I < 1400, and as for larger £ valnes. In order to compnte the signal to noise, 
(4.3), one reqnires the bispectrnm for all £i e [2,2000]. We achieve this by compnting a grid 
of valnes and fitting to a basis in order to smoothly interpolate to arbitrary confignrations. 
The method nsed to achieve this fitting is described in Appendix A. 


4.2 Intrinsic and Doppler Bispectra 

The Fonrier space bispectrnm is defined as 

(,5(ki)5(k2)<5(k3)) = {27rfS^^\ki + ka + kYB{ki,k 2 ,k 3 ), (4.10) 

where the Dirac delta condition (dne to statistical homogeneity) imposes the triangle con¬ 
dition, snch that the bispectrnm may be expressed nsing only the wavennmbers, k^. The 
bispectrnm for the intrinsic (i) and Doppler (d) terms are denoted for a e {i,d}. Using 
eqnation (2.15) and re-writing the Dirac delta fnnction as (27r)^(5^^(k) = J d^xexp(ik • x), 
the CMB bispectrnm for the intrinsic and Doppler terms is inferred to be given by 

^ f f dkidk2dk3{kik2k3YB^‘"\ki,k2,k3)lY^i^j^f {kix dec) jtiikix) , 

(4.11) 

where a e {i,d} labels the intrinsic or Doppler terms, and we introdnce the symbols e 

and to simplify notation. The line of sight integral, J dx, is carried ont between last 

scattering and today, while we choose the npper limit of the integration with respect to the 
wavennmbers ki to ensnre convergence (more explicitly we choose femax = 0.05 Mpc”^). Com- 
pntation of the CMB bispectrnm is, in general, a nnmerically prohibitive task taking 0(10^®) 
compntations for Planck resolntion of ^max = 2000, nnless the Fonrier space bispectrnm B^°‘^ 
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can be expressed in a separable form. In Appendix A, we describe a mechanism for perform¬ 
ing a separable decomposition, given a grid of valnes. The mechanism may also be nsed to 
decompose the ISW bispectrnm. This allows for an efficient calcnlation of the CMB bispec- 
trnm in terms of a basis Qn of separable terms with coefficients for a e {i,d, ISW}. The 
signal to noise may then be compnted readily nsing eqnation (A. 11) as with the correlations 
between the different terms. 

In this section we ontline the calcnlation of the intrinsic and Doppler Fonrier bispectra. 
We neglect the mixed terms which are expected to be negligible following onr considerations 
of the power spectrnm. Following the Green’s fnnction techniqne in §2.3, and mirroring the 
steps nsed in §3.2 for the power spectrnm, we calcnlate the bispectra according to 


B^^\kuk2,k3) = 


(0+(ki,r/i)0+(k2,r/2)0+(k3,??3)) 


V 


(4.12) 


The three point nneqnal time correlator is given in the Ganssian model by [15] 
(0+(ki,r/i)0+(k2,r?2)0+(k3,?73)) _ /3o 1 - i^l2/kU \ 

V " ee^klk^-Kl^''" I n/ 6 j"" \ n/6 ) 

xexp — (^ 2^12 + (^1 - « 12 Mi)^ 13 ) + 5 perms of ki , 

V 6 / 

(4.13) 

where /3o = 8//^(29'7r/3)(F®/?) and r]ij = rji - rjj. In the above formnla we set ^ - 

assnming, without loss of generality, that rji < r/ 2,3 -j with the factors on the second line 
acconnting for nneqnal times, rjij t 0. 

5 Results and Discussion 

We evalnate the GMB bispectrnm dne to cosmic strings indnced by the intrinsic, Doppler 
and ISW terms separately. We neglect cross-terms in onr analysis. Onr goal is to assess the 
comparative improvement achievable in signal to noise from an ISW only bispectrnm - as 
considered by Planck - once last scattering effects are acconnted for. 

The compntation is carried ont nsing the formnlae given in §4. We assnme the cosmo¬ 
logical parameters listed nnder eqnation (2.16). Onr estimates for the signal to noise, (4.3) 
reqnire the calcnlation of the angnlar power spectrnm (7}°*, inclnding the effects of noise and 
beam effects. We firstly generate the inflationary GMB power spectrnm^, nsing GAMB 
[30] and inclnde Planck-iealistic beam, b^, and noise, W by setting = C'^^ + A'^/6^ - nsing 
the prescription ontlined in [31] - where 

bl = exp(-£(£ + l)(T) 3 eam)i '^beam = / i = f^pix^pixj ^pix = — • (5-1) 

VO In 2 iVpix 

9b denotes the resolntion of the beam, Apix labels the nnmber of pixels of area Dpix, reqnired 
to cover the sky, with giving the variance per pixel. In the case of a Planck-realistic 
experiment we nse the parameters 9b = 7.1', Upix = 2.2 x 10“® and Dpix = 0.0349. 

^One should note that our calculation of the cosmic variance represented in the denominator of (4.3) 
neglects the cosmic string contribution, which is known to contribute at a level of less than 3% [3]. 
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The bispectra b^^J^^^,B'^{ki,k 2 ,k 3 ),B‘^{ki,k 2 ,ks ) are evaluated on a grid of 200^ values 
for ii e [2,2000] and ki e [lO”"^,0.5]Mpc~^ (with the latter logarithmically spaced). We 
perform the separable decomposition of the shapes, as described in Appendix A, and, - in the 
case of the intrinsic and Doppler terms, - perform the integrations necessary to obtain the 
corresponding CMB shape. We finally obtain the coefficients of the basis expansion of 
the CMB shapes 


(a) _,(a) {2h + iy/% 2 i 2 + iy/\2i3 + iy/^ 


/TTotTftatTltot 

W2 W3 


Y.ay^Qn{h,i2,h), 


(5.2) 


where a e {i,(i,ISW}, and Qn represents the 3D basis of separable shapes (see Appendix A). 
In this paper we use 160 partial waves, such that we achieve an accuracy of greater than 99% 
accuracy in our decomposition. 

In order to visualise our results, we plot a series of isosurfaces of the quantity 
- normalising with respect to the constant Sachs-Wolfe bispectrum (so called as it is the CMB 
bispectrum resulting from a constant primordial shape S{ki,k 2 ,k 3 ) = (A:iA: 2 /c 3 )^i?(A:i, A; 2 ,/cs) 
with the transfer function approximated by its form in the Sachs-Wolfe limit) - whose form 
is given by [19] 




1 


(2£i + l)(2£2 + l)(24 + l) 


1 


1 


+ ^2 + 1^3 + 3 ii + £2 + is 


(5.3) 


In Figure 2 we plot the individual bispectra for the ISW, intrinsic and Doppler, noting the 
difference in structure between the plots which indicates a relatively low correlation, with 
the Doppler bispectrum of opposite sign to the other contributions. In Figure 3 we plot 
for reference the local model bispectrum - observing the ‘blobby’ features in the plot due to 
the presence of acoustic peaks and the characteristic peaking of the local shape along the 
edges of the ‘tetrapyd’, i.e. at the corners of the triangles on slices of fixed scale .^1 + 1'2 + £ 3 - 
As may be observed, the cosmic string signatures are easily distinguishable from primordial 
shapes such as the local model; due to the active seeding of structure, the acoustic pattern 
is largely washed out - even for the intrinsic and Doppler terms from recombination. The 
total bispectrum is also plotted in Figure 3. The shape is largely negative, as for the ISW 
only bispectrum. Due to the elimination of an £~^ scaling using the constant Sachs-Wolfe 
bispectrum, we observe that the bispectrum is largely of a constant type shape (for £ < 1500) 
with moderate squeezing in the edges (indicated by the deeper blue shade) - for larger £ 
values the bispectrum would be expected to exhibit an additional drop off. The shape 
is compared to the local bispectrum to emphasise the distinct signature - with little evidence 
of acoustic peaks, as expected. 

In order to assess the relative magnitudes, we calculate the signal to noise (4.3), obtaining 
the 1 a bounds 


[A(G/i/10-®)3]* = 0.16, 

(5.4) 

[A{Gfi/10-y]'^ = 0.33, 

(5.5) 

[AiGfi/lQ-^yf^ = 0.15. 

(5.6) 


Note that the bounds are of the cube power of since the bispectrum is proportional 
to {G^y. By contrast the la bound used by the Planck team in obtaining their signal is 
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Figure 2: Clockwise from top left we plot the ISW, intrinsic and Doppler bispectra (nor¬ 
malised by the constant Sachs Wolfe bispectrum). The three shapes are clearly distinguishable 
with the Doppler term of opposite sign to the ISW and intrinsic bispectra. Note that the 
colour bars are not normalised equally. 



[A(G/u/10“®)^]pianck = 0.2. Therefore, the semi-analytic treatment used in this paper for the 
ISW signal is consistent with the detailed - and computationally highly involved - treatment 
used in [3]. 

Using the correlation measure (A. 12) we hnd the following correlations between the three 
contributions 

C(Intrinsic, Doppler) =-0.39; C(Intrinsic, ISW) = 0.03; C(Doppler, ISW) = 0.03, (5.7) 

demonstrating the low correlation between the last scattering surface contributions and the 
ISW bispectrum. Combining the ISW, intrinsic and Doppler bispectra we obtain the bound 

[A(GAt/10'®)^] = 0.11. (5.8) 

Thus we conclude that the power spectrum - for which the 95% bound G/i/10“® < 0.15 is 
achieved [3] - remains a superior probe to the bispectrum when last scattering surface effects 
are accounted for in addition to the ISW only effect. 
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Figure 3: (left) Bispectrum of the local model bispectrum (normalised by the constant Sachs 
Wolfe bispectrum). The acoustic peaks are evident along with the characteristic peaking of 
the shape along the edges of the shape, (right) We plot the sum of the ISW, intrinsic and 
Doppler bispectra (normalised by the constant SW bispectrum). The shape is seen to show 
a mild peaking along the edges, but is largely constant (i.e. proportional to the constant 
SW bispectrum) and negative. The comparison to the local model emphasises the distinct 
signature of the cosmic string bispectrum, with the comparative absence of acoustic peaks 
most obvious. 


One may consider whether the bispectrum may become a superior statistic for detection 
at higher resolution, where the GKS effect dominates the cosmic string signal. In order to 
estimate if this is the case the signal to noise due to the power spectrum can be compared to 
that of the bispectrum. More precisely once compares 


/S\2 ^ (2£+l) Cj ^ 

Uj(2) t 2 (Cft)2’ Uj(3)^ 3 ^ 


(5.9) 


where we simplify the signal to noise (4.3) using + 1)^/^ fa 0.117 for valid con¬ 

figurations^. Assuming for large £ - in the noiseless limit - that the power spectrum, (7^°* is 
dominated by the cosmic string power spectrum, Cg, and using the fact that at high i (> 3000) 
the cosmic string power spectrum and bispectrum are given by the ISW only contribution 
(such that l^Ci = constant and = constant), we find that 




0 . 16 £max 


(^^^ue)‘e=3ooQ 

(^^Q)f=3000 


(5.10) 


and deduce from (3.6) and (4.7) that the signal to noise of the power spectrum dominates that 
of the bispectrum for large values of £max where the power spectrum itself may be dominated 

^One may be concerned that - in this regime - the cosmic variance noise for the power spectrum (and 
the bispectrum) may have contributions from the trispectrum. In [13, 14] it was shown that - while the 
contribution the trispectrum contribution may be of comparable magnitude to the power spectrum (squared) 
contribution - it does not dominate. Therefore, the assumption of Gaussian errors is unlikely to lead to 
erroneous conclusions and is adopted here for the purposes of estimation. 
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by the cosmic string contribution. This analysis has assumed a noise-free signal, while in 
reality the spectrum will be contaminated by point source signals, and Sunyaev-Zeldov’ich 
effects, which would greatly complicate the analysis. Unless such contributions correlate 
strongly with the cosmic string signal the conclusions presented above are not expected to 
change. 

6 Conclusions 

In this paper we have included the last scattering surface contribution in a calculation of the 
CMB bispectrum induced by cosmic strings, for the hrst time. To do this we have used a 
simplihed treatment. First, we modelled the source 3-point energy-momentum correlations 
using a Gaussian model of the string network, which has also been used to calculate the ISW 
bispectrum [11, 14]. The 3-point correlators were integrated with the Green’s function for 
a coupled matter-radiation fluid in the tight-binding approximation. For simplicity we used 
a flat sky approximation, such that the bispectrum obtained is not expected to be accurate 
below I < 200. 

We checked our approach by comparing the angular power spectrum it produces with 
that of a full Boltzmann code treatment, using the Unconnected Segment Model to model 
the string network [21]. We have verihed that on scales I > 200, the qualitative and quantita¬ 
tive features of the power spectrum are reproduced, giving conhdence to the efficacy of our 
approach. 

The main goal of the paper was to assess the signal to noise of the last scattering 
surface bispectrum with last scattering surface effects included, compared to the ISW-only 
contribution. The ISW-only bispectrum produced by our approach compares well to the 
ISW bispectrum computed from maps of the microwave anisotropy produced by the Gott- 
Kaiser-Stebbins effect around a numerical simulation of Nambu-Goto strings [3, 32, 33]. Our 
calculation of the bispectrum contributions uses the separable basis decomposition technique 
- outlined in Appendix A. This allows us to very efficiently calculate the signal to noise of the 
various terms. 

We obtain the the Icr bound A(G/x/10'®)^ = 0.11, with the ISW only bound 
A(G/r/10“®)jg-\Y = 0.15 consistent with that obtained with detailed and highly involved sim¬ 
ulations performed by the Planck team [3]. With cosmic strings constrained using the GMB 
power spectrum to G/i/10“® < 0.15 at 95% conhdence, we conclude that the bispectrum re¬ 
mains an inferior probe to the power spectrum with which to search for cosmic strings in 
Planck data, even with the last scattering surface contributions included. 

Finally, we examine whether the bispectrum might take over as a better statistic at 
higher resolution, where the sharp edges in the maps of the GKS effect are better resolved. 
By comparing the signal-to-noise of an ideal cosmic string GMB power spectrum with the 
string bispectrum, we show that constraints from the power spectrum are always superior to 
those from the bispectrum. 
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A Separable Decomposition 

In the case of the intrinsic and Doppler terms, caicniation of the CMB bispectrnm nsing 
eqnation (4.11) is an almost intractable task, in general, nnless the Fonrier bispectrnm is 
separable, i.e. nnless the bispectrnm, B(ki,k 2 ,k 3 ) can be written as a snm of terms of the 
form f{ki)g{k 2 )h{k^), for arbitrary fnnctions of one variable f,g,h. In addition for the case 
of the ISW bispectrnm, while one may compnte the angnlar bispectrnm, ^ series of 

conhgnrations, caicniation for all £i e [ 2 , 2000 ] is also prohibitive. 

Fitting the ISW bispectrnm from a grid of valnes to an appropriate (separable) basis 
allows one to infer the bispectrnm at arbitrary conhgnrations (alternatively one may ntilise an 
interpolation rontine). In [16, 19, 20] a scheme was developed to allow for an arbitrary (bnt 
snfhciently smooth) bispectrnm to be decomposed into a snm of separable basis fnnctions, 
satisfying onr reqnirements. We briehy snmmarise the approach nsed here (for more details 
refer also to [34-36]). 

Fourier Bispectrum Decompositions — The Fonrier shape is dehned as S{ki,k 2 ,k 2 ,) = 
{kik 2 k‘i)‘^B{ki,k 2 ,k‘i), while the Dirac delta condition ensnres that the bispectrnm is dehned 
on triangles for hxed valnes of ki + k 2 + k^. Utilising appropriately dehned one dimensional 
polynomials of order n, qnikjk^sj^) (where fcmax is the maximnm wavennmber considered) 
- one may nse Legendre modes, or, as done in this paper, those dehned in [19]. The three 
dimensional basis, Qn, is given by triples of these polynomials, labelled by an index n = 
{ni, 712 ,^ 3 } which are partially ordered. In particnlar, we dehne 

Qn{ki,k 2 ,k 3 ) = I;{qni{kilkinax)qn2{k2lkmsx)qn3{k3lkma,x) + 5permntationsof {nj} . (A.l) 
6 

One may dehne an inner prodnct. This measnre is especially nsefnl since it allows for the 
dehnition of a correlation fnnction. An obvions dehnition is to integrate over all allowed 
conhgnrations. The triangle condition sets this region to be 

Vk = {ki,k2, ks'-Y^ki^ 2max{ki),ki ^ /cmax} • (A.2) 

i 

In the case of the intrinsic and Doppler bispectra we choose /cmax = 0.5Mpc“^ to ensnre 
convergence of the integration to calcnlate the CMB bispectrnm. The inner prodnct is given 
by 

(/ 5 ) = / dkidk2dk3f{ki,k2,k3)g{ki,k2,k3), (A.3) 

with the correlation between shapes / and g given by Corr(f,g) = {fg)l\/{ff){gg)- Setting 
din = (SQn) and 7 ^^ = {QnQm), One may approximate the shape, S, as 

S{ki,k2,k3) = Y,anQn{ki,k2,k3 ), (A.4) 

n 

where an = T,m Clearly the ht to the shape nsing the basis expansion will be improved 

with more modes. In this paper we nse 160 modes, for which the ht to the shapes is fonnd to 
be greater than 99.5% nsing the correlation measnre dehned above. 
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CMB Bispectrum Decomposition — The Fourier decomposition allows the projection to the 
CMB bispectrum, (4.11), for the intrinsic (i) and Doppler [d) terms to be performed in an 
efficient manner using 

f + 5 permutations of {n,} j , (A.5) 


where a e {i, d} and the functions are given by 

= - f dkk'^qn{klkjnax)j^^'^\kxdec)je{kx) 

TT J 


(A. 6 ) 


where we recall that e and Writing equation (A.5) in the form = 

En otnQn we note that computation of the signal to noise (4.3) is still a numerically challenging 
task since Qn is not an a priori separable function of the (.i once the integration over the line 
of sight, X, is performed. While one may utilise this form, and calculate the signal to noise 
from each combination of the line of sight times (for each bispectrum factor), leaving the line 
of sight integrations to the end, a more efficient scheme was developed in [16, 35] for carrying 
out this computation for general bispectra. This involves carrying out a late-time (multipole 
space) decomposition. As for the Fourier case, one defines an inner product, given in this 
case by 


ifg)) E ff-MidiMi 1)1/3(2^2 + l)^/^(2^3 + 1)1/3 


(A.7) 


where hg^g^g^ is given under equation (4.3). The denominator in the final factor in the inner 
product such that the weight in the inner product measure is almost constant for all config¬ 
urations satisfying the triangle condition. Inspired by the signal to noise measure, (4.3), and 
the inner product, (A.7), one defines the CMB shape, sg-^g^g^, in the form 




(24 + l)i/®(24 + l)i/®(24 + 1)1/® 

^Cg,Cg,Cg, 


bixi^h ■ 


(A. 8 ) 


To simplify notation we set + l)i/®/\/C^], such that we may write 

~ Sn ctn'^^0n(-^i)4)4)- Given the inner product one may replicate the steps in the 
previous section. One wishes to decompose the CMB shape of the intrinsic and Doppler 
bispectra in the form 


= (A.9) 

n 

such that the shape is in explicitly separable form. Performing the CMB inner product on 
both sides with the basis Qm, results in the expression 

E ^nhnm = E , (ATO) 


where we use the notation = {{QnQm)) and = {{Qri'^ Qm))- Therefore, the coefficients 
are easily found by applying the inverse 7 ,^ to equation (A. 10). 


19 





This form is useful since the signal to noise, (4.3), is given succinctly by 



/sky 


E 


-(a) E 


-(a) 

a; ' 


nm'^m 


(A.ll) 


For shapes, 5*-“^ and with coefficients and an \ respectively, one may also define the 
correlation measure, C, 


C{s 


(a) gib)^ _ 


—(fl.) ^ —(b') 




\/y a^^\/y 

V Linm^n inm^rn y L, 


nm Inm^rn ^ " " 


.(A.12) 


Note on the ISW Bispectrum Decomposition — Unlike the last scattering surface terms, the 
ISW bispectrum is expressed directly in multipole space using the Limber approximations 
of the flat-sky formula (2.12). Having calculated it on a grid of li values we scale it to 
obtain the bispectrum shape, using equation (A. 8). Given that the measure of the 

multipole space inner product is constant for almost all configurations, we instead use the more 
computationally efficient inner product, (A. 3), (where ki ii) to perform the decomposition, 
labelling the coefficients of the expansion to be with = T.na^n^Qn{h,i2 ,h)- 

The signal to noise is then obtained using equation (A.ll) as for the intrinsic and Doppler 
terms. 
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